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Abstract 

In this paper we study analytical solutions for fermion localization in Randall-Sundrum (RS) 
models. We show that there exist special couplings between scalar fields and fermions giving us 
discrete massive localizable modes. Besides this we obtain resonances in some models by analytical 
methods. 
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INTRODUCTION 


The study of gravity resonances in Randall-Sundrum models BE] is of great importance 
because, since stable states for most fields do not exist, it points to the possibility of finding 
unstable massive modes in the membrane. Much of these studies used smooth versions of the 
background with the RS warp factor recovered for large The emergence of resonances 

is an interesting phenomenological possibility that has been approached mainly numerically. 
However, some examples have been studied where an integrable solution was found for 
bosonic helds mm In order to obtain these solutions a specihc form of the background 
was chosen such that a simple potential was obtained for an associated Schrodinger equation. 
That leads to the question: what kind of background allows analytical solutions? The 
answer to this question was given in ra where the authors found an associated equation 
leading to the correct background. That was possible by the introduction of an arbitrary 
function describing the region near the (thick) brane. This arbitrary function makes the 
connection with the effective potential of the Schrbdinger like equation. With this, the 
backgrounds used in the literature mmm are just particular solutions of this equation. 
Besides, many others solutions were found. In particular, the fermion case was studied in 
ref. [H] but the solutions were found only numerically. It happens that the fermion case 
introduces additional complications because of the coupling with the scalar held. Anyway, 
it is possible to circumvent this and generalize the idea used for the gravity case: here we 
hnd a few examples of analytical solutions. Therefore, the complete structure of resonances 
associated with each of these new solutions are also obtained. The way chosen to describe 
the background of the thick brane gives us a better understanding of its structure and can 
lead to new possibilities. 


This letter is organized as follows: in the hrst section we give a brief review of spin-1/2 
fermion background. In the second section we describe fermion in a prohle with piecewise 
warp factor. The third section is devoted to study fermion in a background with smooth 
warp factor. Finally, we present the conclusions and perspectives. 
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THE FERMION SETUP 


We regard here the problem of studying fermions in a Randall-Sundrum like scenario. 
Usually it is considered the coupling between fermions and the gravitational helds alongside 
a coupling between fermions and a scalar held 0, the stuff the membrane is made of. Indeed, 
the membrane is regarded as a domain wall embedded in D = 5. Studies of fermion local¬ 
ization in this context reveals a very natural mechanism of producing chiral helds in D = A 
[lH6l [T 5 H 23 ] . The new ingredient here is that, if the warp factor is dehned in a piecewise way, 
so must be the scalar held 0. The usual coupling is the Yukawa one that will be generalized 
to allow a more complete analysis. 

First, we give a brief review of the fermion background. Let us consider a single scalar 
held thick brane (j) with a potential V{(j)) in (4-|-l)-dimensional curved space-time with 
cosmological constant A. The action describing this model is given by 


^ 


1 (ij _ 2A) - - V{4,) 


( 1 ) 


where R is the scalar curvature, g = det[( 7 Mv], M,N = 0,1, 2, 3,4, and nl = SnG^ with 
G 5 the hve-dimensional Newtonian gravitational constant. We will use the metric in the 
conformal coordinates and the line-element is assumed to be 

= QMNdx^dx^ = e^^^^\'q^ydx^dx'' + dz'^), (2) 

where is the warp factor. The metric tensor presents signature in the form 

(—, -|-, -|-, -|-, -|-), and stands for the coordinate of the extra spatial dimension. The scalar 
held is considered to be a function of 2 ; only: 0 = (j){z). From the metric ansatz ([^ the 
Einstein tensor is given by 

Gfj,u = (3y4^^ -|- 3A")rif^ig , G 44 = 6A'^ , (3) 

Gifi = 0 , (4) 


where the prime denotes derivative with respect to the coordinate. 

The Euler-Lagrange equations for the action Q with the metric ansatz ([^ after some 
manipulations are given by: 


= 3(^/2 _ ^ 

u(0) = + 3A")e-2^ - - , 

^^ = (3AV' + 0")e"'^- 

d(j) 


(5) 

( 6 ) 

(7) 
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If the Z 2 symmetry is imposed, namely A{z) = A{—z), we can see from Eqs. 0 and 0 
that in the limit \z\ —)■ 00 with (f)'{z) —)■ 0 and ld(0) —)■ 0 we have A{z) —)■ — ln(fco| 2 ;| + (3) and 
A = — Gfcg. In this case the bulk is asymptotically AdS. It is worthwhile to mention that the 
condition limi^i^oo —)• 0 can be satisfied by kink solution (lim 2 _,.±oo 0 (^) = ±co, cq > 0 ). 
We shall restrict the bulk asymptotically AdS and kink solutions (with (j){z) = —(j){—z)) 
through this paper. We also make K 5 = 1 for simplicity. 

In this background, the Dirac fermionic action for a spin 1/2 held coupled to the scalar 
held 0 is 


^ 1/2 = J d^x^g , ( 8 ) 

where Dm = Om + xjm being the spin connection. The T^ = are the Dirac 

matrices in the hve dimensional curved space-time and are the vielbeins: = 

gMN_ 

The Dirac equation in hve-dimensional space-time can be derived from the action0 

[Ydf, + 7 ® {d^ + 2d^A) - T]e^F{(j))] T(a;, z) = 0 , (9) 

where is the four-dimensional Dirac operator on the brane. 

Due to the presence of 7 ® in the Eq. g, is convenient to make the chiral decomposition 
of the fermions in Kaluza-Klein (KK) modes [H]: 

^{x,z) = [lpLn{x)'ll:'l{z) +'ll:Rn{x)'lpl{z)] , 

n 

( 10 ) 

where ' 0 Ln(a;) = — 7 ^' 0 Ln(a^) and ' 0 Rn(a;) = l^i^Kn^x) are the left-handed and right-handed 
components of the Dirac fermion helds on the brane, respectively. They satisfy the 
four-dimensional massive Dirac equations in the form of 7^9^t/Ln(a;) = and 

■y^d^'ip^nix) = mn'ipLnix). The KK modes 'ip±{z) obeys the coupled equations 

^ ( 11 ) 

From these we obtain the Schrodinger like equation for each mode 

+ U±{z)'il:{z)± = m^'il:±{z), ( 12 ) 
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where the effective potential U{z) depends on the warp factor as bellow 

U±(z) = (ijf 


( 13 ) 


The study of the massive modes must be done more carefully. This happens because 


when the potential in (12) has the asymptotic behaviour lim| 2 |_>.oo = 0, '4>±{z) are 

linear combination of plane waves and therefore not normalizable. 

For massive modes, like in non-relativistic quantum mechanics, we will use the trans¬ 
mission coefficient to study possible resonances [9]. Usually, the transmission coefficient 
represents the probability flux of the transmitted wave relative to the incident wave. In our 
case, the transmission coefficient represents the probability of a fermionic massive mode to 
pass through the brane. The idea of the existence of a resonant mode is that for a given 
mass the transmission coefficient has a peak at this mass value. In our study of resonances 


we limit the solutions of (12) that have asymptotic planes wave form. 


THE PIECEWISE CASE 


= 


In this section the warp factor A{z) will be dehned in two regions, \z\ < d and \z\ > d\ 

-ln(/(2)), \z\<d] 

(14) 

^-\Yi[{ko{\z\+/3)], \z\>d, 

where d, beta are positive constants and f{z) is called associated function with f{z)>0 in 
\z\ < d, fco = -y/—A/ 6. This is enough to establish how the membrane curves the spacetime 
with the bulk being asymptotically AdS[Ill 112] ■ 

In order to get the desired smooth versions, the functions A, A' and A” must be continuous 
to give 0 and cf)' continuous at \z\ = d. That restricts the form of f{z) in the membrane. 
For the gravitational held the choice of m was A{z) = I lncos(\/l^| 2 :|) giving the effective 
potential U = —Vq for \z\ < d. That allows a detailed analytical study of the model. 
Recently, this result was found as a particular one belonging to a wider class of new exact 
solutions for the gravitational held 0. The purpose here is to present a new way to 
generalize this idea in order to obtain analytical solutions for the fermion held as well. 

We focus our attention in right fermions. The analysis is similar for left fermions. In 


order to implement the boundary conditions, we restrict to even functions f{z) in (14), 
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guaranteeing that the boundary conditions are satished in both edges of the brane, z = d 
and 2 ; = —d. The condition A(0) = 0 implies that /(O) = 1. Since f{z) is an even function 
we also have ^'(0) = 0. Using the equations ([^ and (|^ we obtain 




q/"(z) \ \ ^ j 

3^, 1^1 <d, 


= 


0 , 


and 


(16) 


\z\ > d, 


v{<p{z)) = -r{z)f{z)-6r{zr+6ki 

From Eq. ( [IT| ) we can see that the massless right modes obey the condition 

fo{z) 

With this condition, the potential of the schrodinger-like equation becomes 


U(z) = 


Mz)' 


(16) 


(17) 


(18) 


If we choose fo(z) such that (12) is analytically solved we can obtain analytically the reso¬ 
nances. 

In order to implement the boundary conditions, we are restricted to even functions /o(^) = 


fo{—z) with U{z) = U{—z) in (18). This guarantees that the boundary condition in z = d 
is satished also in z = —d. Our method consist in hnding U{z) and F{(j){z)) that have 
analytical solution for Eq. ( 12 ). As we can see from 0 . F{(l){z)) is an odd function in z. 


From Eq. (15), (j){z) = cq for z > d. We also have two regions for fo{z): 

go{z), 1^1 < d, 

fo{z) = < 


(19) 


{\A+PV^ \A>d^ 

where 7 = riF{co)/ko and go{z) being an even function. In order to obtain a localizable zero 


mode, 7 < —1/2. Now the potential in (18) becomes 


9q b) 
90 (d ’ 


\z\ < d, 


U{z) = 


( 20 ) 


T'U-b 71 > d 

(|2|+/3)2> FI 7 a- 
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The solution of (12) for \z\ > d is a linear combination of Bessel functions of order I 7 I + 1/2. 
With the above considerations and by imposing continuity of A[z), A'[z), A''{z), /o(^) and 
/g(z) at z = id we obtain 


/(d) — ko{d i (3), 
fid) = fco, 

r (d) = 0 , 

9o{d) = (d i (3y, 
9o{d) = 7(di/3)'^"\ 


with the conditions /(d) > 0 . 


The solution for right fermion of (12) will be given by 


= < 


i CmEmy) 1 ^ i d, 

Am9m{,y i \z\ ^ d, 

DraEmi^Eji -2- i d, 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


(26) 


where Em{z) = Em{z) = with u = m{\z\+/3), v = (I 7 I il/ 2 ) and 

9 m{z){hm{z)) are the even(odd) functions. The functions Em{z) and Em{z) have asymptotic 
behaviour In order to analyze the resonances, we fixed the coefficient of the incoming 
wave to one. 

Using the fact that gm{z), Ejn{z), Em{z) are even and hm{z) are odd and the continuity 
of the ipm{z), y'raiz) at \z\ = d we finally obtain the transmission coefficient 


T(m) = 


m‘^\W{g^, h^){dy 


(27) 


\W{E^,gJ{d)W{E^,hJ{dW^ 

where W{f,g){z) = f{z)g'{z) — f'{z)g{z) is the Wronskian. As applications of the method, 


consider the scalar field in \z\ < d given by (/(^) = aA/3sin(z). From Eq. (15) we obtain 


- / {z) + a cos iz)f{z) = 0 


(28) 


The even solution of (28) is given by the Mathieu cosine functions [25]: 




( 29 ) 
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/(z) = C(-l/8,1/6, z) 



FIG. 1. The Mathieu cosine function C'(—1/8,1/16,2;). 


with f{z) real and /(O) = 1. From Eqs. (21), (22) and (23), we obtain d = 'k ko = 0.65447, 
13 = 1.02593 for a = 1/2. The scalar potential can be written as 

^(0) = ^(1 - ^,arcsin(20/\/3)^ 

-6C"^ ^,arcsin(20/A/3),^ + 6kl (30) 

where C\n, m, x) is the derivative of C(n, m, x) with respect do x. The function C{a, b, x) is 
obtained from the real part of MathieuC{a, b, x) divided by the real part of MathieuC{a, b, 0), 
where MathieuC{a, b, x) is dehned in the Mathematica software. For completeness, we show 
the plot of f{z) = C'(—1/8,1/16, z) in Fig. 


For go{z) = AQCOs{^/Voz),Vo > 0 in (19) we obtain U{z) = — Vq, \z\ < ti/2. From Eqs. 


(24) and (25) we have 


7 = -(2 + 

Tlo — (^ + /;^)^/cos(,/V^-). 


(31) 

(32) 


Since 7 < — 1 / 2 , we have some restrictions for the values of Vq. We also have that 
cos(A/ld)|) 7 ^ 0 and sin(v^f) 7 ^ 0. This implies that ^/Vo cannot be an integer. For 
Vo = 1/4 we have 7 = —1.29836 and Aq = 0.409675. We show in Fig. the exponential of 
two times the warp factor and the massless right fermion wave function for Vq = 1/4. 


Now using the equation (17), we show below the scalar fermion coupling 


(33) 


The solution of (12) for \z\ < d is a linear combination of gm{z) = cos(\/m^ + Vqz) and 
hm{z) = sin(\/m^ + V 02 ;). Now using the Eq. (27), we show in Fig. the transmission 
coefficient for the right fermion. As physical result, we have a massive resonant mode near 
m = 0.9. 





















FIG. 2. (a) The exponential of two times the warp factor, (b) The zero mode of right fermion for 
a = 1/2 and Vq = 1/4. 



FIG. 3. The transmission coefficient of right fermion for a = 1/2 and Vq = 1/4. 

THE SMOOTH CASE 

In this section we consider a smooth warp factor given by 

A{z) = ln(/co^^ + 1), -6/co = A. (34) 

This type of warp factor was used in [26] to obtain resonances for Kalb-Ramond helds. From 
Eqs. ([^ and (|^ we obtain 


4>{z) = '\/3arctan(fco2^), 
y{4>) = ^cos2(—). 


(35) 

(36) 


We consider for now, the Sine-Gordon-like fermion coupling: F{(f)) = sin(0/\/3). With 
this we obtain the potential 


U{z) 


koV {-klz"^ + ko'r]z‘^ + 1 ) 

{koz^ + 1 )^ 


(37) 
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Making the substitution 'il){z) = <h(^)(l + and x 

confluent Heun equation 123 


—k^z"^ in (12) we arrive at the 


x{x - l)w"{x) + ((^ + ^)^ “ \w(^) + “ x)w{x) = 0, 


(38) 


where *h(^) = w{x). The right massless mode /o(^) = <^ 0(1 + is normalized for 


rj/ko < —1/2. The solution of (12) for massive right fermion is 


i,{z) = (AHii-klz^) + BH^(-klz^)) (1 + 


(39) 


where Hi{x), H 2 {x) are the confluent Heun functions[2^ and A,B are constants. In this 
case, we show the transmission coefficient in Figj^ for /cq = 1 and r] = —1. With these 
parameters we do not have massive resonant modes. 



FIG. 4. The transmission coefficient of right fermion for Sine-Gordon coupling, with kQ = l,rj = —1. 


Discrete massive modes 


Looking at the expression of the potential, given by Eq. (13), we can study what type of 


F(0) coupling gives discrete massive modes. Asymptotically, for any warp factor we must 
have A{z) ~ 1/z. In order to get bounded states, U{z) must go to a constant or to infinity. 
In the case of U{z) asymptotically going to a constant, we must have that F{(j)) ~ z. For 


the warp factor given by Eq. (34) and F{(p) = tan(0/\/3) = k^z, the right potential i 


IS 


U{z) = 


koV 


klvj'^z'^ 


+ 


{1 + klz‘^YF l + fc2^2- 


(40) 


In this case, the right massless mode is given by fo{z) = Cqs’^o jg normalized 

for 1 ] < 0. For the potential in Figj^ we obtain discrete massive modes m = 2.16, 2.67, 2.87 
for fco = 1 and r] = —3. 
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FIG. 5. The right potential profile for F{4i) = tan((/)/\/3), with fco = = —3. 

It can be shown that is possible to obtain analytically discrete localizable modes. As 

mentioned before, the massless right mode obey T]e^F{(j)) = /o//o, then the right and left 
potential reads 

t/+ = (41) 

= (42) 

From the arbitrariness of we can choose fo{z) in order to obtain normalizable modes 

for "0+ and 'ip-. Let ns consider the case where f/+ — t/_ = 7 , with 7 constant. Then we 
obtain the solntion 

h(z) = Cei‘\ (43) 

U±(z) = ± (44) 

(45) 

Since fo{z) is the zero mode for right fermion, the localization occnrs only if 7 < 0. In this 
case, taking 7 = — 2 n;^, we obtain the schrodinder-like eqnation 

-ip'!^{z) + T oj‘^)'ip{z)± = m^'i(j±{z). (46) 


This is the exactly the Schrodinger eqnation for the one dimensional harmonic oscillator. 
The massive modes are localizable only if m = a/^cj, for right fermions and m = y/2,n + 2 CJ, 
for left fermions, n = 0,1, 2 ■ ■ ■. This gives the conpling 

sin(0/\/3) 


riF{(j)) = 


h 


'0 COS' 




(47) 


This type of conpling was stndied in another context in |28] for a symmetric dS branes. 
This conpling provides analytical treatment and, in this particnlar case, gives a remarkable 
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characteristic: localizable massive modes. It is then important to discuss this coupling in 
order to hnd some physical interpretation, if it exists in fact. Despite this, massive fermionic 
modes can be found in the context of inflating baby-skyrmion brane models [211130] and in 
other more usual scenarios EH. 

CONCLUSIONS AND PERSPECTIVES 

Here we have shown that it is possible to generalize the formalism to obtain analytical 
solutions to models of localizations of fermions in RS scenarios. This work is a generalization 
of previous results only applied to bosonic helds, namely, the gravitational, vector and scalar 
helds. The main idea is to make use of a warp factor dehned in a piecewise way and solve 
an associated equation that describes its prohle centered at the origin. With the analytical 
solution at hand it is quite easy to obtain and study the resonances in the models. We 
analyzed the conditions of fermion-scalar coupling which give localizable massive modes. 
For a very specihc case we found an analytical coupling between fermions and the scalar 
held which gives localizable discrete massive modes. For a very specihc case we found an 
analytical coupling between fermions and the scalar held which gives localizable discrete 
massive modes. These results are very interesting since the mass spectrum of particles is 
discrete. We believe that, despite the fact we are studying fermions, this formalism could 
open a new way to approach open problems as analytical study of non-abelian gauge held 
in branes|32]. More interesting yet is the discussion of models containing another kind of 
spinors, a result which can be applied to supersymmetric models. All of these questions are 
left for future works. 
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